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^ ■ Abstract 

\^ • We determine the u and d quark contributions to the proton magnetic form factor at finite 

momentum transfer by applying chiral corrections to quenched lattice data. Heavy baryon chiral 

perturbation theory is applied at next to leading order in the quenched, and full QCD cases for the 

"^ I valence sector using finite range regularization. Under the assumption of charge symmetry these 

I ' values can be combined with the experimental values of the proton and neutron magnetic form 

^^- factors to deduce a relatively accurate value for the strange magnetic form factor at Q^ = 0.23 



j^. GeV2, namely G%^ = -0.034 ± 0.021 fiN- 
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Strange quark contributions to the properties of the nucleon have attracted a lot of 
interest since the originally puzzling EMC results concerning the proton spin [l|. While 
that motivation has faded [21, y], it remains a central issue in QCD, especially with respect 
to lattice QCD, where such terms necessarily involve so-called "disconnected graphs", i.e., 
quark loops which are connected only by gluon lines to the valence quarks. Despite enormous 
effort y^, the direct lattice calculations of these contributions have so far been unable to 
produce a result which differs statistically from zero. On the other hand, by using the 
constraints of charge symmetry, which is expected to be accurate at the 1% level or better p, 
|6|, one can write relations (c.f. Eqs. flT7|) and flTSl) . below) for the disconnected contributions 
to physical form factors [7|] in terms of valence quantities, which can be accurately calculated 
in lattice QCD, and the experimentally determined form factors. In the case of the strange 
magnetic moment and charge radius of the proton, this approach has succeeded admirably [8|, 
|9|. Here we apply the technique to the strange magnetic form factor at Q^ = 0.23 GeV^. 

Parity-violating electron scattering (PVES) has proven to be a valuable tool for experi- 
mentally determining the strange quark contribution to the electromagnetic form factors of 
the proton. Under the assumption of charge symmetry, one can deduce the strange electric 
or magnetic form factor {G% j^{Q'^) ) from measurements of the corresponding proton and 
neutron electromagnetic form factors and the neutral-weak vector form factor of the pro- 
ton, through its contribution to PVES. While PVES measurements are very challenging. 



a number of groups have succeeded, starting Avith SAMPLE at Bates [10(] and then A4 at 
Mainz [lli] and GO [12I and HAPPEX [H, lli, [l5| at Jefferson Lab. A global analysis of 
all this data has given very precise values for the strange quark contribution to the proton 
magnetic moment, as well as its charge radius [16], which are consistent with the theoretical 
calculations mentioned above. The motivation for our current work is the knowledge that 
in the near future we expect new measurements from A4 and GO at Q^ = 0.23 GeV^. 

In addition to the extensive experimental activity, a variety of theoretical models have 
been applied to the calculation of the strange nucleon form factors. These approaches in- 
clude the QCD equalities supplemented with constituent quark model assumptions [iTJ, 
heavy baryon chiral perturbation theory [l8|, [l9|, dispersive approaches |20L |2l|, |22|, vec- 
tor dominance model (VDM}_J23|], VDM with a kaon cloud contribution 2J], the Skyrme 
model [2 
models l30| 




the NJL model [26|], the chiral soliton |27|, l28 



chiral bag [2 

a two-component model with a meson cloud [33| , etc 



and chiral quark 
These theoretical 



predictions vary quite widely. For example, the predicted strange magnetic moment varies 
from relatively large and negative, —0.75 ± 0.30 [17,] to sizeably positive, -1-0.37 [29l] . 

As well as the above model calculations, there have been some lattice simulations of 
the strange magnetic moment, with early lattice simulations giving a relatively large neg- 
ative value [tI, |3J, |35J. In 2003, Lewis et al. j^ used low order, quenched chiral perturba- 
tion theory, together with the lattice QCD simulation to calculate the strange form factors 
from lattice data. The magnetic form factor which they obtained at Q^ = 0.1 GeV^ was 
-1-0.05 ± 0.06. Recently, by combining the constraints of charge symmetry with new chiral 
extrapolation techniques and low mass, quenched lattice-QCD simulations of the individual 
quark contributions to the magnetic moments of the nucleon octet, a precise, non-zero value, 
G'^(O) = -0.046 ± 0.019, was obtained |8i]. 

In this paper, we present the lattice prediction for the strange magnetic form factor at 
Q^ = 0.23 GeV^. We first extrapolate the u and d quark contributions to t^he proton magnetic 
form factor in quenched, heavy baryon chiral perturbation theory [i 



37j . The quenched 



lattice data from the CSSM Lattice Collaboration is used and finite range regularization 



(FRR) is applied in the extrapolation, because of its improved convergence behavior at 
intermediate and large quark mass [38|, |39|, |40|, |4l|, |42|. In the following we briefly introduce 



the chiral Lagrangian which is used in the extrapolation. The formal calculation of the 
magnetic form factor is then explained, followed by the numerical results. 

There are many papers which deal with heavy baryon chiral perturbation theory. For de- 
tails see, for example, Refs. |43|, |4J, |45| . For completeness, we briefly introduce the formalism 



here. In the heavy baryon chiral perturbation theory, the lowest order chiral Lagrangian for 
the baryon-meson interaction, which will be used in the calculation of the electromagnetic 
magnetic form factors, including the octet and decuplet baryons, is expressed as 

£, = tTTB,{vV)B, + 2DTTB,Si^{A^,B,} + 2FTTB,Si:[A^,B,] 

-ifi:{v ■ V)T,^ + C{fl^A^B, + B.A^Ti:) , (1) 

where S"^ is the co variant spin-operator, defined as 

S>: = ^7 V^^v . (2) 

Here, v^ is the nucleon four velocity (in the rest frame, we have v^ = (1,0)). D,F and 
C are the coupling constants. The chiral covariant derivative, D^, is written as D^B^ = 
d^B^ + [V^,i?i,]. The pseudoscalar meson octet couples to the baryon field through the 
vector and axial vector combinations 

V, = ^(C9,c^ + C^S.C), A, = 1(C9,C^ - C^^mC), (3) 

where 

C = e'^/^, / = 93 MeV. (4) 

The matrix of pseudoscalar fields, 0, is expressed as 

Xtt° + 4=ri n+ K+ 



Bi, and T^ are the velocity dependent new fields, which are related to the original baryon 
octet and decuplet fields, B and T^, by 

B,{x) = e'"'''^''-''"B{x), (6) 

T^{x) = e*™^^^''^''T^(x). (7) 

In the chiral SU{3) limit, the octet baryons are degenerate. In our calculation we use the 
physical mass splittings for transition meson-baryon loop diagrams. 

In the heavy baryon formalism, the propagators of the octet or decuplet baryon, j, are 
expressed as 

— and — -, (8) 

v-k-A + te v-k-A + te ^^ 

where P'^'^ is v'^v'^ — g'^'^ — (4/3)5*^5'^ and A = rrij — itln is the mass difference between the 
baryon j and nucleon. The propagator of meson j (j = vr, K^ rf) is the usual free propagator: 

(9) 
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FIG. 1: Leading and next-to-leading order diagrams for the proton magnetic form factors. The 
last diagram, c, need only be included in the quenched case. 



In the heavy baryon formalism, the electro-magnetic form factors are defined as: 

< B{p')\J,\Bip) >= u{p') \v,Ge{Q') + ''"'''^''"^^''^ Gm{Q')]u{p), (10) 

where J^ is the charge current, q = p' — p and Q^ = —q^. In this paper, we focus on the 
magnetic form factors in each quark sector, aiming to extract the strange quark contribution. 

With the Lagrangian given earlier, the leading and next-to-leading order diagrams for the 
magnetic form factor are shown in Fig. 1. In full QCD, the first diagram, a, is the leading 
diagram, while diagram b gives the next-to-leading order non-analytic term, because of the 
mass difference between octet and decuplet baryons. The last, or so-called double hair-pin, 
diagram need be considered only for the quenched case, where the t] is degenerate with the 
pion. 

The contribution to the magnetic form factor of Fig. la is expressed as 

Gy<2=) . z^ /a. »_(^ + 7/2>^(5-7/2) ^ 

8^'/' J uj{k + 'q/2yuj{ k - -q/2y 



^j{^) = V "^i + A; 2 is the energy of the meson j. We regulate the loop integral using 

finite range regularisation, with u{k) the ultra-violet regulator. Both the pion and koan are 
included in the calculation. In full QCD, the coefficients are obtained from the Lagrangian. 



In the quenched case the coefficients are obtained as in Refs. [36|, |4J, |46 



The contribution to the magnetic form factor of Fig. lb can be written as 
Gb (n2) = d^M! f ^^^Ku(k +^/2)u(k -^/2){uj(k +uj{-^/2)+u(k -^/2)) 



87r3/2 J A 

(12) 
where 

A = uj(k +^ /2)uj{l^ -^ /2){uj{l^ +^ /2)+^){uj(k -^ /2)+^){uj(k +^ /2)+uj(k -^ /2)). 

(13) 
In the preceding equations, /3*(z = a, 6) depends on the quark type, meson loop type and 
whether the calculation involves quenched or full QCD in the calculation. 






FIG. 2: Feynman diagrams at the quark level, which are included in Fig. la for the proton magnetic 
form factor. 



In the quenched case, the additional double hair-pin term from the rj is expressed as 



G\AQ' 



(3F - D) 



'MlGu^Q' 



2887r3/2 



rf^A;- 



"M 



^{kY 



(14) 



where Mq is the double hair-pin interaction strength. We note that the integral of Eq. (14) 
gives rise to a logarithmic divergence in the chiral limit. As a result we estimate the con- 
tribution of this graph using the renormalized value of Gm(Q^) obtained from the lattice 
simulation results at finite quark-mass values. Of course, in full QCD no such term need be 
included. 

In the above formulas, the coefficients in quenched, valence and full QCD can be obtained 
with the same method as in Ref. [46|. For example, the diagram Fig. la is shown in detail 
in terms of the underlying quark lines in Fig. 2. In quenched QCD, the diagram with a 
quark loop has no contribution. In the case of valence quark sector, as well as the quenched 
diagram, the diagram with quark loop can also have contribution if the external photon field 
couples to the valence quark. In full QCD, both the valence and sea quark (loop) couple to 
the photon field. For the pion loop, in full QCD, Fig. 2a and Fig. 2c give contributions, while 
in the quenched case. Fig. 2a and Fig. 2b give contributions. The coefficients for Fig. 2c and 
Fig. 2i are the same as Fig. 2e, which is known from the Lagrangian, since QCD is flavor 
blind. For the same reason, the coefficients for Fig. 2d and Fig.2h are the same as Fig. 2f. 



TABLE I: Coefficients, /J"^, for quarks in quenched, valence and full QCD for Fig. la. The left 
three columns are for an intermediate n meson and the right three columns are for an intermediate 
K meson. 



quark 


u 


d 


s 


u 


d 


s 


Quench 
Valence 

Full QCD 


-ID- 

-4F2 - Id' 

-{D + Ff 


4n2 
-^D- + ADF - 2F' 

{D + Ff 










-i(3F + i?)2 AK 
-\{D-Ff SK 
-\{3F + Df AK 
-\{D-Ff SX 



-{D-Ff 

-{D-Ff 






i(3F + D)2 AK 
1{D-Ff SX 



TABLE IL Coefficients, /3 , for quarks in quenched, valence and full QCD for Fig. lb. The left 
three columns are for an intermediate vr meson and the right three columns are for an intermediate 
K meson. 



quark 


u 


d 


S 


u 


d 


S 


Quench 

Valence 

Full QCD 


!^ 

18 

2C2 

9 


18 

2C2 

9 









18 

C2 

18 




h 

9 
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By subtracting the known coefficients from the total coefficients of full QCD, we can get 
the coefficient for each diagram in Fig. 2. The resulting coefficients for each quark for the 
different cases are summarized in Tables I and II. 

As we know, most detailed lattice simulations for the nucleon electromagnetic form factors 
have been computed in the quenched approximation, in which the strange magnetic form 
factor is identically zero. Since the value in full QCD is not large, any direct calculation 
of G\.i will require considerable effort to extract an accurate value. In this paper, we first 
concentrate on computing the contribution of each valence quark to the proton form factor, 
in the physical theory at the physical mass. Then by using charge symmetry and the 
experimental proton and neutron form factors, we are able to extract a precise value of 
strange magnetic form factor using the techniques of Refs. [7|, |8|]. 

The magnetic form factor can be expressed as 



GuiQ^) = flo + a^fnl + a^ml + ^ GV(Q^ 



(15) 



where the parameters Oq, a^ and 04 can be obtained by fitting the quenched lattice data. In 
the numerical calculations, the SU(3) parameters are chosen to be D = 0.76 and F = 0.50 
{qa = D + F = 1.26) and the coupling constant C is —2D. The FRR regulator, or form 
factor, u{k), is taken to be a dipole {u{k) = jYn^rn^, with A = 0.8 GeV), although as 
shown by Young et al. [41| the model dependence associated with other choices is small. 
We use SU(2) chiral symmetry, with only the light quark masses varying and the strange 



quark mass fixed. Thus the i^-ineson mass is related to the pion mass by: 



m% 



-m^ +m^|phy--m^|phy, 



(16) 



which enables a direct relationship between the meson dressings of the magnetic form factor 
and the pion mass. 
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FIG. 3: The contribution of a single u quark, with unit charge, to the proton magnetic form factor 
at Q^ = 0.23 GeV^ versus pion mass. The dotted, dashed and sohd lines denote the quenched 
(finite volume), valence sector and full QCD (infinite volume) results, respectively. 



The contribution of a single u quark with unit charge to the proton magnetic form factor is 
shown in Fig. 3. The dotted, dashed and solid lines are for the quenched, valence sector and 
full QCD results, respectively. The square points with error bars are the quenched lattice 
data obtained by the CSSM Lattice Collaboration [47|. The lattice results were fit with 
finite volume chiral perturbation theory followed by corrections to yield the infinite volume 
results. The FRR quenched chiral perturbation theory describes the lattice data results 
well over the range m^ G 0.1 — 0.7 GeV^. At the physical pion mass, the quenched {^G%j), 
valence i^G^) and full QCD (^Gjtf) values of the magnetic form factor are 1.099 ± 0.165, 
1.221 ± 0.183 and 1.179 ± 0.177, respectively. 

In Fig. 4, we show the contribution of the d quark, with unit charge, to the proton 
magnetic form factor. The three styles of line have the same meaning as in Fig. 3. Again, 
the quenched lattice results are described very well. In contrast with the u quark case, the 
absolute value of the d quark contribution in full QCD is larger than that in the valence 
case. This is consistent with the disconnected contribution and hence the strange quark 
form factor being small and negative. At the physical pion mass, the quenched ('^GJi/f), 
valence (''G'tf) ^^^ ^^11 QCD (jG'l^) values of the d quark contribution are —0.356 ± 0.053, 
-0.383 ± 0.057 and -0.468 ± 0.070, respectively. 
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FIG. 4: The contribution of a d quark, with unit charge, to the proton magnetic form factor at 
Q^ = 0.23 GeV^ versus pion mass. The dotted, dashed and sohd hues denote the quenched (finite 
volume), valence sector and full QCD (infinite volume) results, respectively. 



With the full QCD values of the u and d quark contributions, one can get the strange 
form factor by subtraction them from the proton or neutron magnetic form factor. However, 
because of the small value of G%,j, the error bar obtained in this direct calculation is much 
larger than the central value of G%j. We therefore use the valence contributions, ^Gm ^^"^ 
'"Gf^, which yield a relatively precise value of G%^. 

The proton and neutron magnetic form factors can be written in terms of quark compo- 
nents as [Tl] 
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r<d I if)P 






(17) 
(18) 



M- 



The label / denotes a "loop" or sea quark 
contribution, while the label v means a connected valence quark contribution in full QCD. 
In the equations above, charge symmetry has been used - i.e. the u and d quark contributions 
in the proton are the same as the corresponding d and u quark contributions in the neutron. 
Charge symmetry is known to be accurate at better than 1% where it has been tested, 
primarily in nuclear systems. It has to be assumed in order to extract the strange form 
factors from parity violating electron scattering. Under the assumption of charge symmetry, 
the strange quark contribution in the proton is the same as that in the neutron. 

The contribution from the quark in the loop in Fig. 2 depends only on its mass - i.e. it 
is independent of whether the quark in the loop is labelled m, d or s. The loop contribution 
of each quark can be obtained using Eqs. fITT]) and flT^ with the same coefficients |D^ — 
2DF + 3-F^ and — ^- By calculation of the relevant loops using FRR, we evaluate the ratio 
^i?^ = Gli/^Gij at Q2 = 0.23 GeV^. This yields the value 'i?^ = 0.185 ± 0.038 allowing the 
dipole mass parameter to vary between 0.6 and 1.0 GeV. Then, using Eqs. flT7|) and flTSl) . 



we find 
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(19) 



(20) 



In Ref. [8]], since we were working at Q^ = 0, it was possible to use the measured magnetic 
moments of the nucleon and the hyperons. Since the hyperon magnetic form factors are not 
known at finite Q^, here we must use the extrapolated valence quark contributions (rather 
than ratios) to extract the strange form factor. The experimental values of G^(0.23) and 



G^(0.23) are 



G;.f(0.23) 



0.98 ± 0.01 m and 



GliiO.23) 



0.96 ± 0.01 [49], where Go 



is the dipole function expressed as Gd{Q'^) = 1/(1 + Q^/O.TlGeV^)^. Substituting the 
experimental magnetic moment of the proton (2.793) and neutron (—1.913), we obtain the 
values G^^(0.23) + 2G5^(0.23) = -0.534 ±0.036 and 2G^^(0.23) + G^^(0.23) = 2.075 ±0.041. 
Comparing the latter with twice the value of ^G^ = 1.221 ±0.183, obtained from our chiral 
analysis of the lattice results, it is clear that there is a significant cancellation in Eq. (TT9ll . 
Furthermore, the large value of 2^G^^ means that the corresponding error on G%j{0.23) 
extracted from Eq. (fT9|) will be large. Indeed, we find that Eq. ( TT9|) yields G^j(0.23) = 
—0.083 ± 0.092. (Note that the quoted error bar arises from the errors in the lattice data, 
the experimental magnetic form factors and finally the theoretical uncertainty associated 
with FRR, especially the variation of the mass parameter A.) On the other hand, the 
relatively small value of ^G'lj = —0.383 ± 0.057 means that we obtain a much more accurate 
value of ^^^(0.23) using Eq. ([2U]), namely ^^,^(0.23) = -0.034 ± 0.021. We note that the 
two extracted values of G^j are consistent within their respective error bars and that the 
sign of both, negative, is consistent with the difference between the extrapolations of the 
single quark magnetic moments in the valence and full QCD cases in Figs. 3 and 4. 
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FIG. 5: Theoretical predictions of the strange magnetic form factors. The two lines are for the up 
and low limits of the experimental data with Eq. (j2ip . 



Some theoretical predictions for the strange magnetic form factor are shown in Fig. 5. 
These models give different values of G%i which are all within the current experimental error 
bars. As for the experimental values of G^^, using the same techniques as Ref. [l6|, we find: 



G^(Q2) = 0.044 + 0.93g2 ± v/0.34 - 7.02g2 + 47.8Q4, (21) 

where Q^ is in GeV^. This form, which is the result of a global analysis of all published data 
50|, is valid over the range < Q^ < 0.3 GeV^. 

The issue of the errors in the strange magnetic moment were, of course, a serious issue 



in the earlier paper [8|] , and there and in the companion papers [51|, |52j we explained all of 
the sources of error, including possible charge symmetry violation. The latter led to a much 
smaller contribution to the final error on G%[ than the statistical errors on the lattice QCD 
data. This is also the case here at small but finite Q^ . The dominance piece of the error 
which we quote to G^^(0.23) arises from the errors on the lattice determination of ^G%[ and 
the experimental errors on proton and neutron magnetic form factors, in comparison with 
which the errors expected from all that is known about charge symmetry breaking in nuclear 
physics, namely that it is typically below 1%, really are negligible^. 

To conclude, we have extrapolated the lattice results for the separate valence quark 
contributions to the proton magnetic form factor at Q^ = 0.23 GeV^ in quenched and full 
heavy baryon chiral perturbation theory. The leading and next-to-leading order diagrams 
are considered and all octet and decuplet baryons are included in the intermediate states. 
Finite-range regular isat ion is used in the one loop calculation, both because it improves the 
convergence of the chiral expansion and because it has been shown to permit a connection 
between quenched and dynamical lattice results JSTJ. By using the constraints of charge 
symmetry, we combine the extrapolated d valence quark contribution with the experimental 
proton and neutron magnetic form factors to obtain a surprisingly accurate determination 
of the strange magnetic form factor G^^(0.23) = —0.034 ±0.021. This is the first time it has 
proven possible to extract an accurate value of the strange magnetic form factor at Q^ = 0.23 
GeV^ using lattice QCD results. It will clearly be of considerable interest to compare this 
with the values which will be extracted from the recent A4 and GO measurements at Mainz 
and JLab^. 



^ We note that the size of the potential charge symmetry violation estimated in the calculation of Kubis and 
Lewis [53| is an exception, being an order of magnitude larger than that found in the earlier calculation 
by Miller [5|, [SJ] • These authors used a very large anomalous omega-N coupling, in contrast with what 
we know from NN scattering. In addition, the omega coupling that they use {g^) is much larger than the 
usual one-boson exchange omega-N coupling. We also note that the implications of this work for other 
examples of charge symmetry violation have not yet been worked out. Nevertheless, if we were to use 
their extreme estimate, our result for G|f(0.23) would change from -0.034 ±0.021 to -0.025 ±0.024. The 
difference is very small and in view of the concerns already noted we prefer not to include this estimate 
of the charge symmetry correction in our final result. 

^ The latest measurement of strange quark contribution to the vector form factors was reported after this 



paper was submitted for publication 55|. The new result favors a negative strange magnetic form factor: 
G^^(0.22GeV2) = -0.14 ± 0.11 ± 0.11. 
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